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ABSTRACT 


Dealiased SEASAT SASS vector winds obtained during the GOASEX program 
have been processed to obtain supe rob servat ions centered on a one degree by 
one degree grid. The results provide values for the combined effects of 
mesoscale variability and communication noise on the individual SASS winds. 
Each grid point of the synoptic field provides the mean synoptic east-west 
and north-south wind components plus estimates of the standard deviations of 
these meanfc. These superobservations winds are then processed further to 
obtain synoptic scale vector winds stress fields, the horizontal divergence 
of the wind, the curl of the wind stress and the vertical velocity at 200 m 
above the sea surface, each with appropriate standard deviations for each 
grid point value. The resulting fields appear to be consistant over large 
distances and to agree with, for example, geostationary cloud images 
obtained concurrently. Their quality is far superior to that of analyses 
based on conventional data. 


INTRODUCTION 


.The ability to determine the winds in the planetary boundary layer 
over the oceans of the Earth by a radar scatterometer called the SASS on 
SEASAT has been demonstrated. The objectives of the SASS program were met, 
and perhaps exceeded, as described by Lame and Bom (1982) . A problem was 
identified when it was recognized that errors in the conventional data need 
to be better understood. The sources for the differences between a SASS-1 
winds, a conventionally measured. wind, and a wind derived from a planetary 
boundary layer analysis using conventional data are many and varied. Yet to 
be produced are planetary boundary layer analyses of the winds based on the 
SEASAT SASS-1 data only, augmented by a minimum amount of conventional data 
such as atmospheric pressures at the sea surface and air sea temperature 
differences, The purpose of this investigation is to prepare and analyse 
wind fields from the SASS data. 

The three most important applications of SEASAT- SASS -like data in the 
future will be (1) the correct description of the winds over the entire global 
ocean at an appropriate resolution, (2) the use of these data to produce 
vastly improved initial value updates for computer based synoptic scale 
numerical weather predictions and (3) as shown, by O'Brien, et al. (1982), 
the specification of the wind stress field and the curl of the wind stress at 
the sea surface for oceanographic applications. In this investigation, synoptic 
scale vector wind fields with a known error structure will be produced at a 
one degree resolution from the SEASAT SASS-1 GOASEX data. These synoptic 
sqile wind fields will be used to compute fields for the horizontal divergence 
of the winds in the planetary boundary, the vertical velocity at 200 meters, 

I 

and to determine the error structure of the resulting fields. The vector wind 

stress fields can then be found from the vector wind fields, and the errors in 

these fields can be computed. Finally, the curl of the wind stress can be 

computed, ,with a specified error structure. 

» 

r 

Certain assumptions need to be made concerning the accuracy of the present 
SASS-1 model function and the wind recovery algorithm. Also the drag 
coefficient that relates wind stress to the wind at, ten meters is a matter of 

some uncertainty. If any of these assumptions need modification and updating 
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for future systems, the main results of this study will still be applicable, 

As a -first effort, the atmosphere will bo assumed to be neutrally 
stratified. An application of the Monin-Obukhov theory for non-neutral 

u 

stability as in Large and Pond (1981) would provide an improved wind 

field at 19.5 (or 10) meters for the computation of divergence. The field \ 

for the air-sea .temperature difference would be needed to an accuracy j 

I 

comparable to the neutral stability wind field. f 

* 

f 
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DATA 

The data used in this investigation were provided by the NASA Langley 
Research Center. They consist of the de-aliased SASS-1 vector winds produced 

i 

as the final product of the analysis of the SASS data for the SEASAT orbit 
segments chosen for intensive study during the Gulf of Alaska SEASAT 
Experiment (GOASEX), Preliminary resultss for this experiment are described 
by Jones, et al. (1979). Two workshop reports, Born, et al. (1979) and 
Barrick, et al..(1979) ll plus a 'summary of conventional data analyses, Woiceshyn 
(1979) give additional details. 

'Tie SASS-1 model function was based on the JASIN data with results 
» 

summarized by Jones, et al. (1982), Schroeder, et al. (1982) Moore, et al. 

(1982) and Brown, et al. (1982). As described by Jones, et al. (1982), the 
GOASEX data were reprocessed 'by means of the SASS-1 model function. De-aliasing 
was accomplished by selecting that wind direction closest to the planetary 
boundary layer wind fields obtained from conventional data. The model function 
recovers the effective neutral stability wind at 19.5 meters. 

i 

The winds recovered from the SASS are very densely concentrated over the 
swath scanned by that instrument on SEASAT. Examples can be found for the raw 
data density in Wurtele, et al. (1982). Only pairs of backscatter measurements 
90° apart were used to obtain winds with all combinations of vertically (7) 
and horizontally ( H ) polarized measurements (i.e. V with V , H with H, V with 
H and H with V ) . For this study, no distinction has been made for possible 
effect of polarization. The communication noise errors could be quite 
different for the different polarization combinations. 

t 

.i 

A sample of a data listing for a particular GOASEX pass is shown in 
■ Table la and lb. The data for a north bound pass of revolution 1141 are 
shown. The central part of the full table is repeated in both la and lb. 

In order, the values tabulated in the twenty one columns (as numbered 
with 11, 12 and 13 repeated) are as follows: 

(1) ( Revolution Number, 

(2) the latitude of the mid point of the line connecting the centers 
of the two cells used in the calculation, 

' r (3) the longitude as above, 

I 
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(4) the wind speed in met err/ per second from a boundary layer model, 

(5) the neutral stability wind speed from the SASS, 

(6) the from which wind direction in degree from a boundary layer 
’model clockwise from north, 

(7) the closest alias wind direction from the SASS, 

’’(8)* a quality code for the accuracy of the boundary layer wind, 

(9) a code for the estimate of the precipitation in the area, 

(10) the distance between the two SASS cells, 

(11) the average pf the values in the next two columns, 

(12) the incidence angle in degrees of the forward beam, 

(13) the incidence angle of the aft beam, 

* (14) the pointing direction of the forward beam in degrees clockwise 

from north, 

(15) the pointing direction of the aft beam, 

(16) the backscatter in decibels measured by the forward beam, 

(17) the backscatter for the aft beam, 

(18) ,the noise standard deviation for the forward beam in per cent 
(multiply the antilog of the backscatter by this number divided 
by 100 to get the standard deviation of the measurement), 

(19) the NSD for the aft beam, 

(20) the polarization of the forward beam (0 = H 3 1~V ) and 

(21) the polarization of the aft beam. 

In this investigation, only the elements of the data vector corresponding 
to LAT, LONG, WSP, WDR and FOR AZ were used to derive the results. For some 
portions of some of the SASS swaths, the boundary layer conventional winds 
differed substantially from the SASS winds. The data base for the conventional 
fields will be discussed later. Future scatterometer systems may eliminate 
the need to remove aliases as in Pierson and Salfi (1982), for example. Wind 
fields produced by the methods described herein would then be independent of 
conventional data sources, ex v ,ep't sea surface atmospheric pressure and,, perhaps, 
air temperatures. Sea temperature is remotely sensed and does not vary 
rapidly. 

The GOASEX data set provided by NASA Langley contained the results of the 
analysis of nine passes over the North Pacific. These were portions of 
revolutions 825, 826, 1141, 1183 , 1212 , 1226, 1227, 1298 and 1299 . Extensive 
depictions of conventional data are available for those five orbit segments 
underlined above (Woiceshyn (1979)). 
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THEORETICAL CONSIDERATIONS 


The objectives of this analysis of the SASS data arc to recover the 

synoptic scale wind at a one degree resolution on a spherical coordinate 

grid and’ to derive synoptic scale fields for divergence, vertical velocity, 

the vector wind stress, and the curl of the jrinds stress from these winds. 

For use in a synoptic scales analysis at scsio initial time, the requirement 

is for the winds to b* specified as east-west and north-south components at 

integral intersections of latitude and longitude in tho form of values as 

close as possible to u s C^ 0 > e 0 ) and CA Q , © o ) , where the subscript, s, 

designates an error free synoptic seals measurement with synoptic scale 

gradients Recounted for and their effects removed and with mesoscale fluctua- 

* 

tions and instrument errors reduced, The SASS winds were not measured 
at the location, A o , 0 q , They contain the effects of mesoscale variability, 
and there are errors (sampling variability as an effect of communication 
noise and cell location inaccuracies) , in the measurement of the backscatter 
that in turn result in errors in the calculated winds. 

The locations of the SASS measurements are more or less randomly 

distributed in an area around the desired location, A o ,0 o , but the gradients 

in the wind are systematic and need to be considered. If the SASS-1 wind 

vector recovery algorithm and model function have no systematic bias, the 

mesoscale fluctuations and the effects of communication noise will also be 

random and. have the same probabilistics and statistical properties within an 

area around A ,6 . 

o' o 

The individual SASS winds can be combined in such as way that the effects 
of 'gradients can be greatly reduced as a source of error in finding the wind 
at A o , 0 q , Also the random errors introduced by mesoscale variability and 
communication noise, can be modeled probabilistically, interpreted statistically 
and greatly reduced by means of the application of small sample theory. 

It is- not necessary to separate the combined effects of mesoscale varia- 
bility and communication noise. Their effects can, however, be considered in 


* 

An appendix defines the notation 


I 
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tho interpretation of tho results that are obtained. Often the effects of 
communication noise stand out above the effects of mososcale variability 
especially as a function of aspect angle relative to the pointing angles of 
the ijadar beams of tho SASS. 

f 

Consider a number of measurements, all of the same geophysical quantity 
such as some kind of a Wind represented by u, made in the same area, that 
all ought to bo nearly equal to the some value. These values will not bo 


equal for numerous reasons. Conceptually there is some correct (or true) 
value and the actual values will scatter in a random way about this correct 
value. The measurements can be described as having a probability density 
function (pdf), f(w), with the usual properties of a pdf, namely, 

• f («) > 0 

, 

Cl) 

j 

£(«) dtt a 1 

-CO 


(2) 

and 

J 

f °° 

1 u£(.a)du « u 1 
-00 


(3) 


which is the first moment of the pdf and can be associated with the correct 
value . 

Moments about u ^ as in 

i 

U a U - 


( 03 


yield 

,i * 

j (u - u,) f (w) du ~ 0 

«oo X 


(4) 

and 

J 

/ CO 

C u - u.) 2 £(«) d u- (A u) 2 

-00 


(5) 

Now let 

u = + tA u 


(6) 

or 

i 

‘t = (u - li^/Au 

a 

C7) 




so that f (u) is transformed into f(t). The pdf', f(u) , is transformed into 

* 

f (w^ + tAu)Audt a f (t) dt 




( 8 ) 


I 


* 

and f (t)dt has the properties that 

( w * v 

£ (t)dt *» 1 

</ “00 
w 

t £ (t)dt « 0 

-00 

r M 2 * 

t £ (t)dt a 1 

-00 


ORIGINAL PAGE I" 

OF POOR QUALITY 

( 9 ) 


CIO) 


( 11 ) 


The concept can be generalized if various reasons for the variability 
of the measurements can be identified. Suppose that there are, say, two 
causes for the variability that many differ from one set of measurements to 
another such that 

(Am) 2 * (Am^ 2 + ( Aw 2 ) 2 (12) 

Then the pdf can be generalized to the product of two independent pdfs, 
f(ti)* f(t 2 ), and (7) becomes 

u *» AUj + t 2 Am 2 Cl 3) 


with obvious extentions, if needed, such as, for example, conditions such 
that t^ and t 2 are not independent and covariances are needed, 

Small sample theory can be applied under the assumption that all of 
the measurements are for the purpose of learning more about m, . Suppose 
that M measurements are made as in 

I 

u i ~ U 1 + t li Aw l‘ ! ’ t 2i Au 2 (14) 


where the t^ and t 2 ^ are from (not necessarily the same) pdfs with properties 
defined by (10), (IT) and (12), (£ (t) is the convolution of £(t.) and f(t 2 ))* 


The average value of the u ^ is given by (16) from (11) 


l 


% 


I 
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u * kl u i 


U l * I ^ t li Aw l * M ^ t 2i Aw 2 


° M 1 * M l H A “ - U 1 + M & t (A U l> + CA W 2> 


and the expected value of the average is (16). 


C («) « u x 


The random variable 


*1 ■ k l hi 


has a mean of zero and a second moment of 


M 2 (tp ■ 1/M 0 

so that u c ; lie represented by 

W “ W 1 * *1 CM A«p + t 2 (M Au 2 ) 

= M x + t* |m'* 1 C(Am 1 ) 2 + (Au 2 ) 2 )j J5 (] 

_ 2 

The expected value of w is and the expected value of (w-up is 


S (u - up 2 = M" 1 (Aw 1 ) 2 + (Am 2 ) 2 


From a given sample, it is possible to estimate the mean, u, the standard 
deviation, Aw, and to use the fact that u has a standard deviation with 
reference to the desired true value given by M 1 Aw where M is the sample 
size. 


A number of SASS individual wind values can be combined in such a way 
!• 

as to recover, a single estimate with a greatly reduced variability. When 
these estimates are then combined to form fields, this greatly reduced 
variability, i which is known, can be used to find the variability of various 
derived fields. 


It is well known that the moan of a sample from almost any typical, but 
unknown, pdf will be nearly normally distributed by virtue of the central 
limit theorem so that although the pdf’s of t, t,, and t„ may not be know, 

is “k A 6 

the pdf's of t , and tg will be close to a unit variance zero mean normal 

pdf for most of the values of M that occur in what follows. 

As in any statistical proceedure the moments of a pdf, such as Uy 

always remain unknown. The statistics only provide a way to put bounds on 

the estimate of w, , that is u, that are made narrower because of the sample 

size. In (15), m, M, (A m,)'' + (AmJ = (Am) are all known statistics, where 
2 2 ac 

(Aw^) + (AWg) is an estimate of the variance from the sample. Thus (15) can 

be rewritten as 


- * — J-' 

M. = u - t M 2 Am 


( 21 ) 


where all is known on the RHS except t (the minus sign is not too relevant)]. 
The values of m, when found at a grid of points form a field consisting of 


values of u plus a quantity that provides information on 


the vari ab il ity 


of m. 


As an example, if u is 10.73 m/s, Am is 2 m/s, and M is 25, then about 
two thirds of the time it would be expected that the interval 10.33 to 11.13 
m/s would enclose the true value, Uy If sampling variability quantities 
are kept track of in the finite difference calculations of properties of tne 
wind field, they provide estimates of the sampling variability of these 
properties. 

t 
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The first data processing step was to sort the data shown as an example 

in Tstble. 1 for each revolution into overlapping two degree by two degree sets 

centered on the one degree integer values of latitude and longitude. The 

vector winds found from combining many individual SASS winds around a grid 

point of a model have been called superobservations. A given SASS wind could 

be found in four different sets. As in Figure 1, a given set of SASS winds 

would have latitudes and longitudes that varied from X - 1 to X +1 and 0 -1 

o o o 

to 0 +1, 

o 

A preliminary investigation suggested that the two solution cases and 

the base of the "Y" in three solution cases gave both directions and speeds 

that were systematically different from the other directions and speeds 

within a given two degree square. Since the original data set gives the 

pointing direction of beam 1, say, x o > dealiased winds were checked to see 

if the directions were within 1°. X Plus 89° to 91°. x ulus 179° to 
o o o 0 o - ' A o * 

181 , and X Q plus 269 to 271 . If they were, those particular data vectors 

were removed from the set. In Table la, for example, the data corresponding 

to SASS directions of 288.8, 289.3, 289.6 and 290.3, were deleted. 


When centered on X^ and 0 , the latitudes and longitudes can be. trans- 
formed by subtracting X q and 0 q from each element in the set. The locations 
of each SASS wind will then be defined by values of, AX and A0 that both 
range between plus and minus one. 















? 
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The average latitude and longitude for the data set will be equations 
(22) and (23) for the N + n values in the 2 degree square. 


N+n 


1 MTJi 

>X ,+ AX. “ X + . . " £ AX 
O' 1 o N+n r 


( 22 ) 


B n + A6 t 
o I 


0 + — 
o N+n 


N+n 

I A0, 

1 


J 


(23) 


The values of AX and Ah will usually be close to, but not exactly, zero. 
The location of this point will be in one of the four quadrants of the heavy 
square in Figure 1. By selectively removing n of the data vectors in the set 
based on the location of SX and A0, it is possible to reduce both 5X and A0* 
for the remaining subset below some preassigned minimum as in equations (24) 
and (25) . 


X + AX 
o 

0 + A0 

o 


N 

1_ 

N 


AX. 

l 


N 

V 
6 

1 

N 

l A0 


(24) 

(25) 


This subset of data for each two by two degree square was processed 

further to obtain the results to follow. Also needed for later use will 

2 2 

be the average values of (AX^) , (A0^) and (AX^*A0^). 


WIND VECTORS 

t 

.* 

The N value for the winds in each data set consisted of a speed and a 

I 

(from which) direction (meteorological convention, clockwise from north) . 
The directions were averaged to obtain a value for the average direction as 
in (26) where the summation notation is abbreviated. 


X = | I Xj 


(26) 


Each wind was then resolved into a component in the direction of y 
and a component normal to x as in (27) and (28), 


V p . = |V|. cos (X. - 5) 


(2 7) 
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( 28 ) 



These components were averaged to obtain (29) and (30). 



(29) 

(30) 


In all cases V^‘ was essentially zero. The standard deviations of V p and 
V^j about their sample means were also computed as in (31) and (32) . 


Au = 'SDCVp) = C±- ICVp. - ?p) 2 ) % (31) 

Av = SD(V n ) = (i IV W -.y V* (32) 

k 

The processed data at this point consist of the sample size, N, the 

mean direction, x, the mean component in the mean direction, V p and the 

values of AU and AV as located at a point given by A q + AA (with the values 

of A + AX’ and 0 + A0 located very close to the integer values of latitude 

o o 2 2 

and longitude over the ocean) plus the average values of (AAj) , ( A0j) 

and AA. A0.(i.e. VAR (AA) , VAR (A0) and COV (AAA0) . 


The "from which" meteorological convention can be converted t a "toward 
which" vector by adding 180° (mod 360°) . The new direction, clock angle, is 
needed for use in analyzing the wind field. The east-west, U^, (abbreviated 
as U), and the north-south, Vq, (abbreviated as V) components of the three 
non-zero vectors given by (29), (31), and (32) are next found. 


The mean east-west and north-south components are given by (33) and 

(34). 


0 ss y sin 

(X + 180°) = -Vp 

sin x 

(33) 

V = V^ cos 

(x + 180°) = -V p 

cos x 

(34) 


The components of AU and AV from (31) and (32) in the east~west 
direction according to the convention shown in Figure, 2 are (35) and (36). 


1 


AUj = -Au sin x 
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(35) 


Au, « - Av cos X (35) 

* 

The components in the north-soijth direction are given by (37) and (38) . 
AVj = - Au cos x . (37) 

Av 2 = Av sin x (38) 

The vectors (V p , x) > (AU,. x) and (AV, x + 90°), can b' - ' used to 
two orthogonal vectors, as in (39) and (40). 




Au 


(39) 


V Ns " V N - *2 AV = 0 - h Av t 40 ’ 

and X = X s ‘ The subscript, s, represents the desired synoptic scale value. 

The orthogonal coordinate system has one axis parallel to the direction, x, 
and t^ and t 2 are independent random variables with zero means and unit 
standard deviations as in equation (9) to (11). When resolved into north- 
south and east-west components, east-west variability is correlated with 
north-south variability so as to confine the scatter of points to the form 
of ellipses such as the one shown by the dashed line in Figure 2. The 
variability of 0 and V, where U s and are the desired synoptic scale 
components, analogous to in equation (7), can be represented by (41) 
and (42). The calculated values, U and V, equal the true values lb and V , 
plus the 'effects of gradients and the random errors. 


0(A + AA, 0 + A0) = U (A + AX, 0 + A0 ) 
v o o ' s o *o 

+ h Au i + *2 Au 2 


(41) 


V(X o ^+ AX, 0 Q + A 8) = V s (X q + AX, 0 Q + A0) 

+ t l Av l' + * 2 Av 2 


(42) 


Following these transformations, the data for each integer value of 
latitude and longitude for each point in the SASS swath consists of A , 



e o , AA, A0, 0/ AU r AU 2 , V, LV V AV 2 , VAR(AA) , VAR(A0) , COV(AAA0) and N. 
The following stops process the 0 and V components of the processed SASS 
data to determine the best possible field for the synoptic scale wind. 

'Consider the east-west component of the wind at a point where it was 

* 

estimated by the SASS somewhere within the two degree square. The value 
could be described by equation (43) . 


3 U CX, 6J 

ua„ ej .. u„u„, e o ) -Vi — — u. - y 


s 0' 


SU (A , 0 ) t 

+ ~^ni - y + t t 4U 


where AU is to' be defined later. 
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The average value of the N values of U(A^,0^) is 


0 = i l (Utt^ 8.)) = 0 CX 0 + AA, 0 O + A6) 

3U (A . e ) 8U (A ,9 ) i , 

■ o * -TT- 5 - ^ rr^-® + 


S 0 


a e 


(44) 


Everything except U S (A 0 >9 0 ) is known, or can be estimated, and the two 
other terms on the right hand side are usually small because AX and A0 are 
small.. For points in the interior of the swath, it i s possible to estimate 
a correction from the gradients of U as in (45) where uncorrected values 
near the grid points are used to find approximate gradients, 


0(A o , e o ) = 0(A o + AA, 0 O + m 

f 

I 

- %(U(A rt + 1 + AA, e n + A0) - 0 (A - 1 + AA , 0 + M))AA 

o * o o 0 

- %(U(A 0 + AA, 0 o + 1 + W) - U(A 0 -1 AA, 0 Q ~ 1 + A0))A0 (4o) 

subject to, the conditions that appropriate inequalities hold, namely 
> r 

either U(A b + 1 + AA, 0 Q + AS) > U(A q + AA, 0 Q + M) > 0(A q -1+AA, 0 q +A 0) (46) 
or U(A q + 1 + AA , 0 q +A0) < U(A q +AA, 0 q + A0) < U(A q - 1 + AA,0 q + A0) (47) 
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and either G(A o + 3X,0 +1+Se) > U(A 0 + AA,0 q + M > U(A q + AA, 0 Q -1+ A0) (48) 

or 0(A '+ IX, 0 +1+ A0) < U(X o + Al,O o + A0) < U(X o + AA,0 q - 1 + A0) (49) 

If sets of these inequalities are not satisfied, then U(X qJ 0 o ) is near a 
maximum or a minimum, in one or both directions, (or a saddle point) and the 
correction is simply to set U(A o , 0 q ) equal to U(A o + AA, 0 Q + A0) , 

The expected value of the square of U(A^,0^) - O(A o ,0 o ) is given from 
(43) by (50). 


e C H (U(A., 0 .) - U(A o , 0 o )) 


1 x 
2 


■l- 


i 8 W 0 ^ 


9A 


VAR(AA) + 


90 


VAR (A0) 


3U.(A , 0 ) 8U CA ,6 ) , 2 

+ 2 ~ COV (AA AO) + (AU ) 


3 A 


3 0 


(50) 


In (50), the first three terms represent the effects of gradients of the 
east-west component of the wind. These and similar terms contribute to the 
variability of U and V in (31) and (32) and of AU^, Al^, AV^ and AV ^ in 
(41) and (42) . These effects need to be removed so that the effects of 
sampling variability can be found. 

The contribution of the variability of the values of U(A^,0^) and 

V(A^,0.) from gradients in the synoptic scale wind is the first three 

terms in (45) . For SASS values uniformly scattered over a particular two 

2 2 

degree by' two degree square, the average values of (AA.) and (A0^) will . 
be about one third and the average value of AA^,A0^ will be close to zero. 
SASS winds concentrated in diagonally opposed quadrants would yield a non- 
zero value for this third term. The synoptic scale contribution to the 
variability of the estimate of U s CA Q , 0 Q ) can be estimated by (51) where the 

values corrected by ,(45) have been used. 

» 

r 

VAR (SYNOPTIC) = ~ (0(A + 1, 0) - U(A - 1, 0)) 2 ' VAR AA 
» o o 

+ i (0(X o , e 0 + 1) - O(X o ,0 o - l)) 2 VAR A0 
+ \ CO(X 0 : + i, e 0 ) - 0cx n -i, e o ))(ucx o , 0 o + i) - u(X o ,0 o -D) cov AXAe (si) 


f 


Tho variability of the estimate, OU o ,0 Q ), that is the result of 
mesoscale variability of tho winds from one BASS coll to the othor.and 
of the effects of communication noise is found by means of (52) . 

t 

# 

VAR (MESO PLUS COMMUNICATION NOISE) 

= (Atjp 2 + (AUg) 2 - VAR (SYNOPTIC) 


* CAU* ) 2 . (52) 

where (AU ) is now defined in (43) by equation (52)* Since VAR(SYNOPTIC) 

o o *2 

is an estimate, it may exceed (AU^)" + (AU 2 ) in which case (AU ) is set 
to zero. 'The analyses of superobservations by Pierson (1982, 1983) assumed 
that gradients of the synoptic scale wind could be neglected in the study 
of superobservations. For some parts of the wind fields obtained from the 
SASS, the assumption is not justified, and this correction needs to be made. 

i 

Often it is small. Specific examples will be given later. 

A parallel set of equations for the north-south components of the 
synoptic scale winds can also be obtained. The final results of processing 
the data in this way produces values of u (^ o >0 o )> V(A o >0 o ), AU and AV at 
the integer values of A and 0 in the main part of the SASS swath where AA 
are A0 small. Near the edges of the SASS swath, which will be discussed 
separately, other methods are needed because of the smaller number of values 
in the superobservation, the location of the data fairly far from the desired 
latitude and longitude and the lack of values for all of the quantities 
needed in the above equations. 


The elliptical scatter of the winds that form a superobservation is an 
important aspect of the sampling variability of the measurements. From (52), 
and its equivalent for V and from (35), (36), (37) and (38), reduced values 
of AU^ and AU 2 may prove useful in studying the effects of mesoscale 
variability and communication noise. 


Let 

»» 

K ='^( AU*) 2 + (AV') 2 | 


(AUj ) 2 + (AU 2 ) 2 + 


(AV,)' 


(Av 2 y 


-1 


(53) 




Thou let 
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(54) 


Auj * K* 5 AU X 
AU 2 * AU 2 

avJ * r* av 1 

AV 2 = AV 2 


(55) 

(56) 

(57) 


.This change effectively reduces the elliptical scatter of the original data 
by an amount attributable to the removal of synoptic scale variability over 
the two degree by two degree squaro. 

The fdnal result of the steps taken so far is to make it possible to 
represent the east-west and north-south components of the synoptic scale 
wind at X q ,0 o in forms similar to equations (19) and (21). They are 
equations (58) and (59) and equations (60) and (61). 


UU 0 ,6 o ) - U s U o ,6 o ) + tjN ~ h AU* + t 2 AU 2 (58) 

va 0 ,9 0 ) - V s U o ,0 o ) + A V ; + t 2 N- 51 Avj (59) 

VW = u(x 0 ,e 0 ) - t x n' h Auj - t 2 n'* 5 au 2 ( 60) 

WV = V(X o ,0 o ) - t x A V ; - t 2 N"* 5 AV 2 (61) 


For an analysis of (58) and (59), the same values for t^ and t 2 must 
be used in both equations. This shows that 0 and V are not independent and 
for example that 

1 

e ((t0cx o .B o ) - u s CX o ,e o )) (v(X o ,e o ) - v s (x o .e o )j) (62) 

= N” 1 (AU^ Avj -i- AU' 2 AV 2 ) (63) 

» t 

’ Equations (60) and (61) put all of thekknown statistics on the right 

hand side., It needs to be interpreted with care. There is, conceptually, 

only one correct value for U , and only one, for V . Picking t 1 and t 7 at 

random generates many values of U and V , one of which may be the correct 

s s 

one. It t^‘ and t 2 are constrained to lie on a circle and are normally 


I 
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distributed then 


PC 


U* J. 4 .'" > 

fc . T t ,» S 


'2 






*2 ,2 
*i * *2 


$ 


dt l dt 2 


|R J2n 0 -C 2 /2 w?de 


2tt 

= 1 - o -^ 2 

so that if J? b 2, the ellipse generated in the U V plane will have a 

s s 

probability of 0,865 of enclosing the true value of U g and V g . 

In the analysis of superobservations, each SASS wind in the four 
degree square is given equal weight in contrast to some of the present 
analysis techniques that weight wind and pressure reports from ships as 
a function of how far away they are from the grid point being analysed. 
For cc '.ventional data and conventional analysis procedures, the same 
report will influence a grid point from distances as far away as five 
or ten degrees of latitude or longitude. Other sources of error 
dominate the analysis of conventional data and maximum use must be 
made of each oceanic observation because of the poor spacing and 
sparceness of the data. 

For SASS data the error sources are different and the dominant ones 
are essentially random. Equally weighted observations are the most 
effective way to reduce the variability inherent in the SASS data. 


I 


THE FIELD OF HORIZONTAL DIVERGENCE 

OF P©08 Qii # Adir7 

For constant density at a fixed height above the sea surface, the 
equation of continuity in spherical coordinates is given by equation (64) , 


1? 




. ' 1 3U . 

+ ren 'sx * 


1 8 (cos OVD 

'srri — La 0 


(64) 


The divergence of the horizontal wind at 19,5 meters for a neutral 
atmosphere is given by equation (65) , 


«* 2 "h ° scS* 


os 


+ cos 0 |~ - sin 0 V 


(65) 


A finite difference estimate of the divergence at a one degree resolution 

requires values of U at A +1, 0^ andA *1, 0 and values ofV at A ,0 + 1; 

> o o o o o 

A q , 0 q and A q> 0 Q - 1. IVe neglect also the ollipticity of the Earth, and 
use i? « 2 •lO^Orr)” 1 , The equations are in the form of (60) and (61). 


The finite difference value for the horizontal divergence of the wind 
is given by equation (66) where the various AU^, and so on, are associated 
with the appropriate latitude and longitude. 


(div.wj 
( 2 h s cos 


10 


-6 


(ou o + 1. e 0 ) - tj au; - t 2 iU ; - OCA 0 -1, e 0 ) 

(va 0 . e 0 + n - 1 5 Avj - 1 6 av 2 


1 I . 

+ t 3 Uj + t 4 AU 2 I + cos 0 


- *CV 6 o - 15 + h Av i + h av 2 


- *9 AV 1 - ^0 AV 2 > 


- 0.0349065 sin 0(V(A .0 ) 


( 66 ) 


The expected value of (div 2 W^) simplifies to equation (67) 
* —6 f 

e (div,w.) = — 0(A + 1, 0 ) - 0(A - 1, e ) 

2 h s cos 0_ ( o * o ^ o o ■ 


+ cos e o e?u o , e 0 + i) - vu o , 0 O - 1)) 


0 


- 0.0349065 (sin 0 V(A , 0 3 = div 0 W 

o o' o J 2 h 


( 67 ) 
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The oxpectod value of the variance of the divergence is given by equation (68) 


e (cdiv 2 w h ) s - (div 2 "w h ) i 


« VAR (div^ W^) 


‘ 2.025*10 

$3 *** » » " m il' " ! wn 


-11 


(cos 0 Q ) 


5 - (Auj a 0 + 1. 0 o )i 2 + CAU 2 tt 0 + 1,0 O )) 2 + C&uJ(X o -l,0 o )) 2 


+(au 2 u o - 1, 0 o )) 2 + Coos 0 o ) 2 CCAVj(X 0> 0 o + l )) 2 + (av 2 (X o , G o * l )) 2 


+ C A V ^ o'°o ' 15)2 + “VV V 1))2) 

* 1.2184-10’ 3 (sin 0 o ) 2 CCAV’ (X^)) 2 * CAV 2 CX 0 , e Q )) 2 ) ] .(68) 


The divergence at one of the grid points in the SASS swath where 
sufficient data are available is thus given by equation (69) where t by 
the central limit theorem is approximately a normally distributed random 
variable with a zero mean and a unit variance. 


(div 2 W h > 6 = div 2 W h + t (VAR(div 2 ) )' 1 


(69) 


The wind in the layer of air near the ocean surface can be expressed 
as a function of height, h, given the wind profile for the first few 
hundred meters. Given the friction velocity, (u )» the divergence can be 
expressed as a function of height for neutral stability and integrated 
from the surface upward, This will be done after the fields for the wind 
stress and the curl of thp wind stress are found. 


BOUNDARY LAYER AND WIND STRESS 


The Monin-Obukhov theory and the concept of the drag coefficient given 
by equation (70) for a neutral atmosphere 
*C r 


J D10 


2 ,„2 

u*/U 


10 


where 


u* - t/p - 


i t 

<u w > 


(70) 

(71) 


have been the methods used to try to understand the turbulent boundary layer 
over the ocean. Many different sets of measurements have been made in order 
to try to determine the relationship between wind stress and the wind at 10 
meters. 

i 

It is not relevant to this particular investigation to review the great 
many papers that have been written that describe the results of these many 
investigat.ipns, Reviews or some of these studies are given by Phillips (1977) 
and Neumann and Pierson (1966) . 


More recent results using more modern and more carefully calibrated 
instrumentation that cover a large range of wind speeds are those of 
Davidson, et al. (1981) .Dittmer (1977), Smith (1980) and Large and Pond (1981). 
Low winds are covered by Davidson, et al. (1981), low and moderate winds, by 
Dittmer (1977) and moderate and high winds, by Smith (1980) and Large and 
Pond (1981) . 


» i 

Smith and Large and Pond used methods that measured <u w > directly as 
well as the dissipation in the Kolmogorov range which was than correlated 

i t 

wi ’Jh <u w, > and U^. All of the data obtained to determine the relationship 
between wind stress and U 10 scatter when plotted either as x/p versus U 10 
or as C 10 versus Uj Q . 1 


To investigate the problem in still one more way, an analysis was made 
* 

as a term paper , by Vera, of the data provided by W. G. Large and extracted 

from the publications of Davidson, et al. and Dittmer. For neutral stability, 
* 

the data Were of the form n^, (x/p)^, U^.. , for a data set that ranged from 
light winds* of 1 m/s to winds of 27 m/s. The higher winds have smaller values 


* 

Vera, Emilio E., A Study of Curve Fitting Procedure for the Wind Versus 
Wind Stress Relationship. Term paper for an Oceanography Course, The City 
College of New York. 


for n. (whore n. represented the number of individual runs in a restricted 
i i 2 

range of wind speeds that were averaged to produce the u* and U 1Q values) 

and weight the fit less strongly than moderate and low winds. 


; fhe*need is for the most accurate possible prediction of the wind stress 

i 

given the wind at ten meters. Instead of first finding C 1Q , it seems that it 
would be more direct to minimize the error in predicting the wind stress from 
the wind speed from the available data as in equation (72) . 


N 

Q = E(n.(t/p) i - A U 


lOi " B U 10i 


C u 3 ) 2 
L U 10i J 


(72) 


This is accomplished by finding those values of A, B and C tnat minimize Q 
by requiring that 3Q/3A = 0 3Q/3B = 0 and 9Q/3C = 0 which yields equations 

(73). 


I 'J 101 In. U 10i In 5 U 


l n. U? n . In. U 4 

u i lOi u i 


i lOi 
i u 10i ^ n i U 10i 


£ n i U 10i ^ n i U 10i ^ n i U 10i 


A 


I n. U 10 . (T/p) i 


B 

s 

E n i. U 10i ^ T,,p h 


c 


I n t uJ Qi Ct/p). 

(73) 


A good fit was found for the data that were used in terms of equation (74) . 
u 2 = T/P = 10" 3 (2; 717 U 10 + 0.142 U 2 Q + 0.0761 U 3 Q ) (74) 

, Graphs of x/p versus U^q for the averaged data of Large and as extracted 
from the two other sources are plotted on linear scales in Figure 3 and on 

i 

■" logarithmic scales in Figure 4. Figure 4 shows how the slope of the x/p versus 
Uj 0 relationship varies from one to almost three as the wind speed increases. 

Vera also investigated the possibility that some power law as in equation 
(75) mightffit the data better. 

•» 

x/p = aU P ' (75) 

It was found that no power law did as well as equation (74) . As shown in 


I 


Figure 4 power laws can only be made to fit a portion of the data and fail 

outside of a restricted range. Equation (74) will consequently be used in 

this investigation subj v ct to the comment made in the introduction. Given 

any equation that relates wind stress and the wind at ten meters, the 

analysis that follows can still be applied. 

♦ 

The SASS data are for winds at 19. S meters. To calculate the stress it 
is necessary to know the wind at 10 meters. The wind at 10 meters is 

0(10) = j~~- Cn (10/z o ) (76) 

, since u*(U 10 ) is known as a function of U 1Q , which implies z q , and at IS. 5 
meters it ’is 


U(19.S) 


u*(U 10 ) 

K 


(19,5/z ) 


(77) 


so that 


U *C U 10 ) 

U (19 .5) = U(10) + — — — £n (19.5/10) 


(78) 


When U(10) is varied in convenient increments, a table of U(19.5) versus U(10) 
can be generated. Given a value of U(19.5), interpolation yields a value of 
11(10) to be in the calculation of T/p and if p is known T. 
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THE CALCULATION OF THE STRESS 

p OOR QUALlTy 

Given the individual SASS winds scattered over a SEASAT swath, one 
might, consider taking the basic data, reducing the wind to 10 m., calcu- 
lating the stress and deriving such quantities as the curl of the stress at 
whatever resolution is available. For many reasons, this line of analysis 
will lead to dubious results. 

Aspects of .the difficulty can be illustrated by calculating a super- 
observation wind stress two different ways; The first is by calculating the 
stress at each SASS point averaging the stresses, and analysing the 
statistics, and the second is by using the results of equations (5) to CIO), 
reinterpreting them and calculating the stress and its statistics from the 
vector averaged winds. 


A new coordinate system is needed as a first step. Let two great circles 
intersect at A , 0 with one .in the direction y from (5) and the other at 
right angles to it, and let distance along the first be given by G^ (degrees) 
and along the second by G 2 (degrees) A coordinate transformation would be 
possible as in the divergence analysis as a last step when needed. 

The wind from the SASS at any point in the two degree by two degree 
square would then be represented by (79) and (80). 


V_. = U (A. ,6.) + AU. 
Pi P l l l 


V M . = V M ( A . , e . ) + AV. 
Ni N l i l 


(79) 

(80) 


where in turn 


3 U(A ,0 ) 3U (A , 0 ) 

AG 2i 


+ t. . AU + t_ . AU 
4i m , 2i c 


( 81 ) 


and 


3VCX ,0 1 3V U 0 > 

AV i r ^-5 g: -— AG ii + • as!" - AG 2 i 


+ *31 AV m * ^4i W 3 


( 82 ) 


OF POOR QU&UTY 

In (81) and (82), any residuals from the gradients of the wind from the 
average values of AG^ and AG 2 ^ can be removed for superobservations as 
shown by equation (44) . The terms AU^, AU c> AV m and AV c have a new meaning 
and must be interpreted in terms of the analysis given by Pierson (1982). 

The terms, AU and AV m , represent the contribution to a particular SAGS wind 
of the mesoscale variability of the wind for each SAGS cell. The terms, AU 

r * V 

and AV, , represent the communication noise error in the measurement made by 
GAGG for each pair of cells. They are in the form of standard deviations to 
be multiplied by random numbers - , t^, t^, and t^, from zero mean unit 

variance probability density functions. This random error is a not too well 
known. function of incidence angle, aspect angle, polarization and neutral 
stability \Vind speed. For winds near upwind, downwind and cross wind, the 
communication noise error can be substantial and the mesoscale contribution 
is a strong function of wind speed. If finite differences are used to 
calculate divergences over smaller distances with the original GAGG winds, the 
last two terms in (81) and (82) dominate the calculation of the gradients to 
find differences as opposed to the analysis given previously for a superobser- 
vations. The GAGG G0A6EX data do not provide a way to separate these two 
effects. Even more discouraging results are obtained when these individual 
observations are used to obtain wind stress and wind stress curl. 


From (79) and (80), the magnitude of the wind at is 


|vj = (Cu p O i ,e i ) * AU .) 2 + (AV .) 2 )* 2 


(83) 


and so the magnitude of x/p is 


T/p (A., 0.) * A V. I + B V. 2 + CIV. ' 

r 1 1 ’ x '.i 


(84) 


When expanded to second order in AU^ and AV^, this becomes (85) 
x/pO.,0.) = AUjCA.,©.) + BOJpCAj,©.)) 2 ♦■C(U p (A i ,e.)) 3 


+ (A + 2BU p (A.,e i ). + 3C(U p (A.,0 i )) 4 ) U. 


(B + 3CU p (A i , 0 i )) (AU.) 


+ CB + |cU p (A i ,0 i ))(AV i ) 2 


(85) 
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The squares of AIL and AV^ involve terms of the form, 

t 2 CAU ) 2 + t 2 CAU ) 2 
li ■ m 2x v c 

that do not fluctuate about zero and that do not average to zero. 

Moreover they are amplified at high .winds by the wind speed. Individual 

values of the wind stress calculated from the SASS data thus appear to be 

» 

somewhat untrustworthy. 


If the N winds used to generate the superobservation wind are each 
used separately, N values of T/p can be obtained and these N magnitudes 
of x/p can be averaged. The expected value of this average is given by 
(86) where, U = U(X q , 0 Q ) . 


£ (T/p(X.,e.)) = A0 p + B0p + C0p 

30' 30 

+ Ca+ ‘ 2b0 p + 3 C °pH^- ^ + ^ ag 2 ) 


2 _ , SU p 2 


3U 2 3U 3Up 

+ ( B+3Cu p r VARAG 1 + ^ VARAG 2 + 2 307 3GT C0V AG 1 A( V 

X 4 i X Z, 


3 \ 2 


3V p 2 


3V 2 


3V p 3V p 


+ ^ VAR 4G 1 * ^ VARAG 2 + 2 3G7 C0V 4G 1 AG 2 3 


+ CB+3CUpJ CB+ |cu p ) C-ggi- jg 5 - VARAGj + 

8Up 3V 

+ 3G^ 5 C0V4G 1 4G 2 




VAR AG, + 


3Up 3V p 

3Gi 3G 2 


CAU ) 
m 


CAiy 




cav c v 


( 86 ) 


Although the average can in theory be corrected by various estimates of 
the numerous terms that produce a bias, the errors of various kinds in the 
original SASS data are quite large, Additional terms to third order would 
introduce even larger effects since the cube of the wind speed is involved. 
Variability in direction produces even more intricate results, 

h 

2 

The variance of the data would be found from e (x/p 0^) - x/p (X^,0^)) 
and yields terms involving the squares of AU^, AU c> AV m and AV c multiplied 

t 
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functions of U, A, B and C. This variance would be reduced by N ^ when used 
in calculations involving the wind stress field but proceeding in this way 
leads to intricate correction procedures and the need to correct for large 
effects 

i 

In contrast consider the average wind for the N values around A q , 0 q as 
given by (29) and (30) 'which could be corrected for the effects of a gradient 
by finding (45) and its counterpart for V(A O ,0 Q ). The standard deviations 
caused only by the combined effects of mesoscale effects and communication 
noise could be found from (54) through (57) and returning to vectors parallel 
and normal to the direction y. 

i 

The result would be four numbers at A , 0 q , namely x> and 

T P ■ o’ V + * S N ' !i W mc C87) 





( 88 ) 


where the combined residual effects of mesoscale variability and communication 
• -U 

noise have been reduced by N 2 to describe the sampling variability of the 

mean. Failure to correct for the contribution from the synoptic scale 

gradients and using the values found at A + AA, 0 + AA with standard 

oo 

deviations reduced by N~ 2 would also give much more stable results. 


For a super observation, substitute (87) and (88) into (85) and the 
result is '89) . 

' T/pC V Vs " A0 P - Bfl P * C °P 

I 

+ (A + 2B0 p + 3C0p) t 5 W h AU mc 
+ (B + 3CU p ) t\ N~ 1 (AU inc ) 2 

,r f 

* (8 + fcO p ) t\ N'hiV^) 2 (89) 

Over a one degree field the uncertainty in the value of r/p(A , 0 ) is 

oo 

greatly reduced by the factor N~ 2 and the bias is even more strongly reduced 

i 

by N . Similar effects occur for direction. 
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It appears to suffice to process the data in an even simpler way for the 

center of the swath. The superobservation wind data can be reduced to 10 

2 

meters instead of the individual winds. The value of U* = t/p(A o , 0 q ) can be 
found from 

* 

T/p U 0 , e 0 ; Up) = A Op + B Op + C Up (90) 

The variability of the stress in the direction x can be found from (27), (31) 
and (30) after correction by (54) to (57) and after correction to 10 meters, 

A p (i/p) = T/p (0 p + N“^ AU mc ) - T/p (0 p ) (91) 

i 

The variability normal to x that would give a vector stress with the 
same direction as the vector .sum of (87) and (88) with tg * 0 and t^ = 1 
is given by (92) 


A m (t/p) = t/p(U d ) (- 


N~ 2 A V 


The stress is represented by 


T/PlCVVps = T/pt °P 3 + t 7 V T/p) 


x/p(A ,0 ) M = t 0 A k .(t/p) 
o o Ns 8 N 


with orthogonal components in the direction, x and normal to x« These in 
turn can be resolved into east-west and north-south components and processed 
in, the sariie way as the winds. 


i . * 


TUB CURL OF THE WIND STRESS 
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The curl of the horizontal vector wind stress is found by first multi- 


plying equations (93) and (94) by P(A o »0 o ) which can bo found from the air 


temperature and the atmospheric pressure at the sea surface and then computing 


the north-south and east-west components of the stress plus the appropriate 


Ax's. Values of P(X 0 ,9 q ) vary about 1.25 kg/m 4 . In spherical coordinates, 


31, 


CURL (T^) a ^ 


3 Ceos 0x,* . J 


(A) s' 


a 0 


) 


1 cos 0 . cos e 


a x 


(A)s 


i? 


a a 


a e 


+ sin 0 T U)S ) 


(95) 


where x,„, and x,,. are the north-south and east-west vector components of 

1.0JS V.AJ s 

x . If for example 


T Ce)s ( VV = *<e) c W + S At i * ho At 2 


(96) 


(A >0J = t,,,(A .8 ) + t 17 At, + t, 0 At, 


(A)s 1 o- o J " (A) ^ o' o J T “ll uu 3 T 12 


(97) 


and so on, the finite difference value of (95) is 

-6 


CURL (T h ) s = 


(cos 0 


“1 1 (0)S C V + 1> 0 o 3 " T (0)s ( V l * V 


cos 0 o (t (a)s (A o , G o + 1) - T U)S (A o , e o - 1) 


+ 0.0349065 sin 0 Q x (x)s CA o ,0 Q )) 


(98) 


-2 -1 

‘with dimensions of (newtons m ) m~ . 


Ten different random effects are found in the calculation. The 
expected value of the. curl and the expected value of the variance of the 
cur! are given by (99) and (100) where ten different Ax’s from the five 
different stresses that are used are needed for (100) . 


• i 


* i 
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6 CCURL c?h ” s ■ i'l'ij ' (w** + 1 ’ s ° 3 ■ T «» w ° ■ ij v 

- cos 0 o Ct (x3 (X o ,0 o + 1) - T CX j(A o ,0 o - 1)) 

+ 0,0349065 sin 0 o (X o* ° 0 - , l 
» CURL (T h ) 

VAR CCURLC?.)) J- - f (AT(X + 1, OJJ 2 

h Coos 0 O ) 2 ^ 0 02 


' + catU o - i, e o ) 2 + catcx 0 - i, e 0 ) 2 

+ (cos 0) 2 C(At(A o ,0 o * 1) 3 ) 2 + ) 

+ 1.2184'10 -3 (sin 0 O ) 2 (CAi(A o ,0 o ) x ) 2 + •••)] 

, where •••* the represents the additional terms needed to complete the 
full equation, 

so that 

(CURL T, ) = CURL(x, ) + t (VAR (CURL { t .)))* 1 
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If tho neutral wind profile is assumed to oxtond to the height, h 2 > 
equation (65) can bo rewritten and the result integrated from h^, very near 
the sea surface whore W (I? 0 + h 1 ) is nearly zero at the synoptic scale, to h 2> 
In the limit, h^ can be set to zero and the subcript, 2, can be omitted. 


ft + h 2 


R o + h l 


Ij.cA W)<U? 


ft + h 2 


R o * h l 


7? div 2 W d R 


( 102 ) 


Many of the terms that arise and that enter into the finite difference 
calculation* are negligable, Two turn out to be important and are effectively 
multiplied by terms that are constants as far as the integration is concerned. 
One is of the form given by equation (103) , 

(h 

n 2 

(constp dz B const^ (h 2 - 1^) (103) 

■ h l 

and the other is (104) 

\ h h, 

const 2 en( z/10) dz = constgChg ^ Jq* " h 2 “ h l 071 To" + h l^ C104) 

ft 

By means of equations similar to (27), and following, with the wind 
reduced to 10 meters and with the standard deviations of the means, values 
for Vp, AU and AV can be found at X q , 0 q . 



Define 



n _ itj 

u, - (AV p + BVp + CVp) 

(105) 


Au*p ~ u* (Vp + AU) - (Vp) 

(106) 

and* 

1 * 

=l u„(? p ) • (iv/Vp) 

(107) 

Then 

u* u z/10 AU* In z/10 

Vft =V P + K + (AU + K ) 

(108) 


• v N U) s « t 2 (AV - 


AU* n Ik z/10 


K 


(109) 


The terms analogous to (33) to (40) aro 
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u* Cn (z/10) 

U(z) - -CV p + ) sin x 


’« 0(10) + O) Cn z/10 
uin (z/10) 

V z * -C^p + — ) cos x 


« V(10) + (6^)071 (z/10) 


sin x Au* p Cn (z/10) 
A U lCz ) - AU, ^ 


•* AUj + (fi 2 U x ) An (z/10) 


AU 2 (z) * AU 2 - 


cos x Au^ N %(z/10) 


AV 1 (z) = hV x - 


AU 2 + t<S 2 U 2^ Cz/10) 

cos x Au* p An (z/10) 


= AV 1 + (5^) An (z/10) 


AV,,(z) = AV 0 + 


sin x Au* n Cn (z/10) 


AV 2 + M CZ//10) 


( 110 ) 


(111) 


( 112 ) 


(113) 


(114) 


(115) 


The appropriate terms, all of which can be found for each grid point in 
the field, can be substituted into (71) , After integration of (102) and 
simplification by the results implied by (103) and (104), the final result 
is equation (116), where the subscript on h 2 has been dropped. The All's and 
so on, need to be found at the correct latitudes and longitudes. 
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+ cos 'e o CV 10 CX o> 0 O + 1) + 6 1 VCX o ^ 0 O + 1)C 'fin 


- W Q o -w - «iV(v« 0 " ~ » + 


+ 0.034906S sin 0 (V. n (X ,0 ) + ) 

O 10 0 0 


At h = 200 meters, to ~ - 1 is 2.0 and (117) becomes (118) in ms" 1 . 

i 

*> C200:, s ■ -5STF- °10 C \> + V + 2Cs i UCX o + x > 6 o» + ••• 

O , 

+ cos 0 o (V 1O CX O ,0 O + 1) + 2(5 1 V(X o ,0 o + 1) + ... 


+ 0.0349065 sin 0 q CV 1o CX o , 0 o ) + 26^^, 0 Q )) 


(116) 


4* ♦ • # • 1 


(117) 


As in preceding analyses. 


Lb (200) s = £b (200) + t Alb (200) 


(118) 


which can be evaluated in a way exactly similar to previous results. 

i 

As a note in passing, there is not much difference between (105), (106) 
,and (107) and (90), (91) and (92). The computation of error •; for the 
specific application is made simpler by the choice that was made. From 


.2 . f A. ,2v»% A.. 2-> /, ,2\h 


u* + Au* = (u* + (Au*)) = u*(l + (Au*)/u*) 


= u*(l + (Au*)/2u*) = u* + (Au*)/2u* 


( 119 ) 


it can be seen that two different ways to find Au* yield nearly the same 
quantity. 


INDEPENDENCE AND COVARIANCES 


The reasons for the variability of one SASS pair of backscatter values 
compared to another pair are basically uncorrelated and random. The mean 
wind speed and the mean direction for a superobservation are perturbed by 
two independent and uncorrelated random effects, one parallel to the wind 
and one normal to it. A correlation becomes evident when resolving the super- 
observation into east-west and north-south components. A superobservation has 
been based on data from a 2° by 2° square, and once the effect of synoptic 
scale gradients are removed (which involves SASS values scattered randomly 
over the square from as far away as X Q ± 2 and 0 q ± 2 ) , the estimate of 
the standard deviation of the mean is based on a sample of independent random 
variables. The choice of grid point spacing for superobservations for 
future systems remain a matter to be investigated and is a function of the 
model to be used and of the scatterometer that obtained the data. A 
very coarse resolution model in the horizontal would require a more careful 
treatment of gradients. 

The use of overlapping squares at a one degree resolution introduces 
correlations between the estimates of the winds and the standard deviation 
of the variability simply because the same values are used in different 
estimates. Every other grid point in the east-west or north-south direction 
is independent and those located diagonally (for example, A o ,0 q to a q +1, 0 q +1) 
are weakly correlated. Smoothness in the wind field can be judged by 
inspecting every other point in the grid. 

For the calculation of the divergence, the vertical velocity at h and 

t 

the curl of the wind stress, the analysis becomes more complex. Some 
particular SASS speeds and' directions can make contributions to the means 
plus standard deviations of three of the five terms that enter into the 
calculation. (The additional complexity introduced by (45) is not considered). 

. Equations (69) and (118) are essentially statistics formed from a linear 
»’ 

weighted combination of the original data. Equation (100) is a pseudo- 
linearized plus non linear weighted combination of the original data. 


I 


% 

Correlations do not do much damage in the estimation of means, but they 
can affect estimates of variances and standard deviations. The contribution 
of a particular SASS value to the final estimate of the divergence can be 
traced through the entire analysis, The value for the difference between two 
estimates- of the east-west component so as to find 3U/3A is calculated from 
completely different east-west components of the SASS values and. thus the 
numbers on which the estimate is based are independent, The north-south 
difference so as to find 3V/30 is based on completely different north-south 
components of the SASS values weighed by a cosine and are also found from 
completely different values and thus the estimate is based on independent data, 

For the calculation of the divergence as in (69) and quantities such as 
(101) and (119), the full expression is, however, somewhat more complicated 
and merits a more complete analysis. In Figure 3, the one degree squares that 
contain SASS values that enter into the evaluation of (74) are numbered from 
one to twelve. The values of 0 (A o , 0 q ) and V (A ,0 ) (and so on) are 
partially obtained from SASS values in squares 4, 5, 8, and 9. The wind for 
the point, X q} 0 q , was represented by equations (60) and (61) , Let N equal 
N = + Ng + Ng + Ng, and then the term f t^ AIL can be represented by (120) . 


N 

i 


N, 


N r 


Nr 


“I* l 


A AU. . 
i*4 1,4 


It., 


AU 


1*5 


t i>8 


AU 


1,8 


N c 

- i 

i 


t i,9^ U l,9 


( 120 ) 


where the t's on the RHS are independent and the second subscript represents 
the one degree square containing the SASS value. The expected value of each 
side is zero and the expected value of the square of each side is 


N 4 

£ C I t 
l 


1.4 AU M 


N S 

♦ It ls 


N 


“l.5 + 


8 

l t 


N 9 

1,8 AU 1,8 + t l,9 AU 1,9 3 




( 121 ) 


where the assumption that the error properties are the same in all four 
squares is made. 

All of the contributions from each of the squares numbered in Figure 1 
can be identified by a double subscript notation as above. The analysis can 

• i 

«» 


I 


I I I 

be carried through keeping track of the dependence of Aik and AV- and of ALL 
i — 4. 4 

and AV 2 throughout the finite difference calculation. 


The contribution to the divergence from squares, 4, 5, 8 and 9 come 
from terms involving the gradients of U and V and value of V at A , 0 . 

I | S S S 0 0 

For the AU^ and AV^ terms in (60) and (61) a part of their contribution 
can be represented by ( 122 ) where a is the coefficient of the sine term 
in (67). 


rh 
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sin 
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1,5 


0 
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- cos 
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AV* n 

+ 
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sin 
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AV* n 
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1,9' 
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i..;4 

+ cos 

0 

Av' /i 
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sin 
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Av' . 


0 

1,4 




0 

1,4' 


* h.a N 8 % ( - AU 1,s' - C0S e 0 AV 1,8 * “ Sin e 0 AV 1,8 ] a22) 


The expected value of this contribution to the divergence is zero. The 

2 2 2 2 

expected value of t 1 ^ 9 ^ 4 and t^ g are one. The squares of the 

terms in parentheses give the terms found in (73) subject to the reasonable 

assumption of ( 121 ) and what it implies plus twelve cross product terms. 

The signs of the cross products terms are such that two of four similar tenns 
have positive signs and two have negative signs. The cross correlation 
between 0 and V when estiamted from scattered SASS values will consequently 
cancel out and be close to zero in the estimate of the variability of the 
divergences as in (73). Although the terms in (71) as identified by t^ to 
t ' are not strictly independent, this analysis shows that ( 68 ) and consequently 
(69), is essentially corrept. The same conclusions are also valid for (99), 
(100), (101), (116), (117) and (118). 
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